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Abstract: 

We formulate a consistent expansion of the QCD evolution equations for the twist- 

three quark distributions g2{x,Q^), hL{x,Q^) and e{x,Q^) based on the interpretation 
of the evolution as a three-particle quantum-mechanical problem with hermitian Hamil- 
tonian. Each distribution amplitude can be decomposed in contributions of partonic 
components with DGLAP-type scale dependence. We calculate the corrections to 

the evolution of the dominant component with the lowest anomalous dimension - the 
only one that survives in the large- iVc limit - and observe a good agreement with the 
exact numerical results for Nc = 3. The admixture of operators with higher anoma- 
lous dimensions is shown to be concentrated at a few lowest partonic components and in 
general is rather weak. 
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1 . Twist-three parton distributions in the nucleon are attracting increasing interest 
as unique probes of novel quark-gluon correlations in hadrons with clear experimental 
signature, giving rise to certain asymmetries in experiments with polarized beams and 
targets. Quantitative studies of such asymmetries are becoming possible with the increas- 
ing precision of experimental data at SLAC and RHIC, and can provide for an important 
part of the future spin physics program on high-luminosity accelerators like ELFE, up- 
graded CEBAF etc. With this perspective, a detailed theoretical study of twist-three 
parton distributions in QCD becomes mandatory. 

Altogether, there exist three twist-3 distribution functions - chiral-odd, e(x, Q"^) 
and h^iXjQ'^), and chiral-even, g2{x,Q'^) - each being a function of parton momentum 
fraction x and the energy scale Q^. By virtue of Lorentz invar iance the distributions hi 
and g2 contain contributions of twist-2 structure functions hi and gi, respectively: 

dy ~ dy ^ 

hL{x) = 2x —hi{y) + hL{x), g2{x) = -gi{x) + —giijj) + g2{.x) . (1) 

J X y J X y 

The QCD description of the remaining genuine twist-3 part of these distributions func- 
tions, g2{x) and of e{x), is usually believed to be quite sophisticated. Their 
moments are related through the QCD equations of motion to the matrix elements of 
quark- ant iquark-gluon operators which have a nontrivial scale dependence and mix with 
each other under the renormalization.|T| As a consequence, the QCD evolution equa- 
tions for the functions hiix), g2{x), e{x) cannot be written in a closed form and require 
additional nonperturbative input. An important simplification occurs, however, in the 
large-A'^c limit. It was shown p, ^ that to this accuracy the twist-three distributions 
and e as well as the flavor-nonsinglet contribution to g2 satisfy simple DGLAP-type 
evolution equations (see e.g. 0) 

+ ^'^(1-^) + ^/-^}, (2) 

where cx/ = — 1 , and 1 for / = hL,g2 and e, respectively, and 1/(1 — z)+ = 1/(1 — 2;) — 
5{1 — z) dz' / {1 — z'). This result implies that the inclusive measurements of twist-three 
distributions are complete (to the stated accuracy) in the sense that knowledge of the 
distribution at one value of is enough to predict the distribution at arbitrary Q"^. 
From phenomenological point of view it allows to relate the measurements of different 
experiments to each other and to compare them to the model predictions (including 
lattice calculations) that typically refer to a low scale. 

Although the large- A^^c version of QCD presents a valid theoretical limit, its relevance 
to the actual = 3 world and its accuracy in predicting the scale dependence of twist-3 

t Additional mixture with three-gluon operators is present for flavor-singlet contribution to g2{x,Q^). 
In this letter we consider only the flavor-nonsinglet part. 
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distributions is a priory not clear. This work presents the first attempt to go beyond 
the large- A'^c approximation in a systematic way. In particular, we calculate the 
corrections to the evolution kernels in and show that mixing with quark- antiquark- 
gluon operators remains under control. 

2. From the OPE analysis one finds that the scale dependence of the moments 
of the twist-3 distributions J^^ dxx^'^'^ f{x) is governed by renormalization of the set of 
local composite quark-antiquark-gluon operators {k = 0, N) 

[3^]% = q(D-nf in^\G,,±iG,,^,] (B -nf-'q, 

[TtTn = QiD -nfn^a^P V n^G,, (B -nf-^ , T = {1, ^75} (3) 

for chiral-even (S"^) and chiral-odd distributions {T\ and Tj^^ for e(x) and hii^x), respec- 
tively), see e.g. Here, is a light-like vector and G^u = ^fiupxG^^/2 stands for a 
dual gluon field strength. To leading order, renormalization of Ti and Tj^^ is the same 
and we, therefore, drop the subscript in what follows. Similarly, it is enough to consider 
the operator S~^. 

The operators [T]% (and [>S'^]^) with different k = 0, ...,N and the same number of 
covariant derivatives N > mix with each other under renormalization. The mixing ma- 
trices have been calculated to the leading order (e.g. P) and can be diagonalized numer- 
ically for any given A^. The eigenvectors, then, define the multiplicatively renormalizable 
operators and the eigenvalues give their anomalous dimensions. A disadvantage of this 
(traditional) approach is that the mixing matrix does not have any obvious structure in 
this basis and is not symmetric. As the result, the structure of the spectrum remains ob- 
scure and the eigenvectors are not mutually orthogonal with any simple weight function. 
Relative importance of various contributions is also far from being clear .0 

A systematic l/N^ expansion of the evolution equations is made possible by going 
over to the Hamiltonian formulation developed in in |]^, |^, ^ |TU|, TT|. The renormalization 



group evolution is driven to leading logarithmic accuracy by tree-level counterterms and 
has the conformal symmetry of the QCD Lagrangian. The evolution kernel can, therefore, 
be written in an abstract operator form in terms of Casimir operators of the collinear 
subgroup SL{2, R) of the conformal group. In this way, diagonalization of the mixing 
matrix of the twist-3 operators in (^ can be reformulated as a three-particle quantum 
mechanical problem 

n ^'7V,g(Xi, X2, X3) = £N,g^N,giXi,X2, X3), (4) 

defined by the Hamiltonian 0, ||, [10|, |Tl| 



Ha = N^nf - —n^i\ A = {T,S+}, (5) 



*E.g. the hierarchy of entries in the mixing matrix 0{N), 0(1), 0{l/N), etc. is not preserved in the 
eigenvalues. 
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with 

nf, = v^f ( J12) + ( J23) , Ti^^l = nl'n^i2) + c/i? ( J23) + c/,? ( J13) . (6) 

Here, the notation was introduced for two-particle quark-quark and quark-gluon kernels 
V;(°)(J) = ^|;(J + 3/2) + ^|;(J - 3/2) - 2^|;(1) - 3/4 , (7) 

U^^\J) = i,{J + 1/2) + J - 1/2) - 2V'(1) - 3/4 , 

«^ (J-3/2)(J- l/2)(J+l/2) ' «^ 2(J-l/2)' 

y«( J) = V'( J) - V'(l) - 3/4 , J) = \ [i;{J - 1) + V'( J + 1)] - V(l) - 3/4 . 

Here and below %l}{x) = d\nr{x)/dx stands for the Euler ^-function; the subscripts '1,2,3' 
refer to antiquark, gluon and quark fields, respectively. The operators Jik are defined as 
follows 

Jik {Jik — 1) = Lik — {J'i + J'kT 1 (8) 

where Lj arc the generators of the SL[2, R) group and L^^ are the corresponding two- 

— * 

particle Casimir operators. The generators Lj can be realized as differential operators 
acting on coordinates Xi of the wave function ^(xi, X2, x^) 

L_^i = Xidl. + 2jid^. , L+,i = -Xi , Lo,i = Xid^. + ji , (9) 

where the conformal spin j is equal to ji — js = 1 for the (anti)quark and j2 = 3/2 
for the gluon field, respectively. The eigenvalues of the operator Jj^ define the possible 
values of the conformal spin in the two-parton channel and are given by ji + jfe + n with 
n being nonnegative integer. 

The Hamiltonians defined above are manifestly SL{2, R) invariant: 

[n,L^]^[n,L^]^[L'',L^]^0, (10) 

where La {a = 0, -|-, — ) are the total three-particle SL{2) generators 

La — La,l + La,2 + -^a,3 , — Lq[Lq — 1) + L+L_ . (11) 

One can, therefore, diagonalize the three operators L^, Lq and Ti. simultaneously. The 
conditions 

Lo-^N,q{Xl-,X2.,X^) = {N +l/2)-^N,q{Xl-,X2.,X^) , L_-^ ^ ,q{Xl-, X2-, X^) = Q (12) 

define ^'jv,g(xi, ^2, ^3) to be a homogeneous polynomial of degree N that does not contain 
factors of {xi -\- X2-\- x^) and therefore does not vanish as Xi — 0. 
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Solving the Schrodinger equation one constructs the basis of multiphcatively renor- 
mahzable local operators, O^^g. Omitting the Lorentz structure, the correspondence is, 
schematically 

ON,g = ^NADq,Dg,Dg)qGq (13) 

with covariant derivatives Dg, Dg, Dg acting on the antiquark, gluon and quark, respec- 
tively. The corresponding eigenvalues provide the anomalous dimensions £N,q- 

where h = llNc/3 — 2nf/3. Here, parameter q enumerates different eigenstates of the 
mixing matrix. 

Note that the conformal operators defined in ([13| ) diagonalize the full mixing matrix, 
including mixing with the operators containing total derivatives. Since at the end only 
forward matrix elements of these operators enter the parton distributions, taking into 
account this additional mixing may be seen as an unnecessary complication. The ad- 
vantage for such formulation is, however, that the Hamiltonian becomes hermitian with 
respect to the scalar product 

(^i|^2) = / Vxxixlx3%{xi)'^2{x^) (15) 
Jo 

with Vx = dxidx2dx36{xi+X2+X3 — l). As well known in quantum mechanics, hermiticity 
of the Hamiltonian implies that all its eigenvalues (anomalous dimensions) S]\f^q are real|^, 
and the eigenfunctions are mutually orthogonal: 

{'^N,q\'^N',q') = SNN'Sgg'\\'^N,qf, \\'^N,qf= Vx Xixlx^]'^ N,q{Xi)\'^ ■ (16) 

Jo 

Completeness and orthogonality of the eigenstates corresponding to multiphcatively renor- 
malizable operators will be of crucial importance for this work. In particular, using these 
properties one can expand the moments of the twist-3 parton distributions over the con- 
tributions of multiphcatively renormalizable operators as 

f'dxx^^' fix, Q^) = Y: ^-^^{{OnAQ'))) ■ (17) 

j'l g V^^A 

Here, {{ON,q)) are reduced (scalar and dimensionless) forward matrix elements of the 
corresponding local operators over the nucleon state and $Ar(xj) are the coefficients in the 

§ Although the conformal symmetry is lost beyond one-loop, it is easy to prove that anomalous 
dimensions of twist-3 operators remain real to all orders in perturbation theory. Indeed, entries in the 
mixing matrix are real and complex eigenvalues may only appear in (complex conjugate) pairs. On the 
other hand, the number of eigenvalues is fixed and in leading order they are non-degenerate (which one 
can establish by inspection). These two conditions are contradictory since a pair of complex conjugate 
eigenvalues (i.e. with the same real parts) cannot be obtained from a discrete non-degenerate spectrum 
at ckg ^ (i.e. with all real parts different) by a perturbative renormalization group flow. 
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re-expansion of twist-3 quark-antiquark operators that define the parton distributions, in 
terms of quark-antiquark-gluon operators. They are given by QCD equations of motion 
and for a;i + 0:2 + ^3 = take the simple form |1TT| 



for f = e ,hL and g2, respectively. Note, however, that in order to form the scalar product 
in (|1^) one needs to know ^N{xi) for Xi + X2 + X3 = 1 that corresponds to taking into 
account contributions to the QCD equations of motion of the operators containing total 
derivatives. The corresponding expressions are given below in (0). Each term in ([T7| ) 
has an autonomous Q^— evolution, Eq. (p!4D, and brings a new nonperturbative parameter 



3. The large-iVc Hamiltonians 7^^'^^ in (^) possess an additional 'hidden' symmetry, 
related to the existence of a nontrivial conserved charge 

2 ^2 2 



[H^\ Qt] — , Qt — {LI2, LI^} - -LI2 - -LI^ , 



[H'§1,Qs+] = , Qs+ = {L\2, -L23} - -L\^ - -^23 , (19) 



where { , } stands for the anticommutator. As the result, eigenstates of Ti^^^ can be 
labeled by quantized eigenvalues g^, (£ = 0, 1, . . . , A^), of the operator Q: 

Q X2, X3) = q ^^xi, X2, X3) . (20) 

Since the number of degrees of freedom (=3) equals in this case to the number of conserved 
charges (L^, Lq, Q), the quantum mechanical system described by the Hamiltonian Ti.^'^^ is 
completely integrable. This allows to use advanced mathematical methods for the analysis 
of the spectrum and, in particular, develop the WKB expansion of the eigenvalues at large 



The remarkable property of the large- Ac evolution kernels P, ^ that is responsible for 
the simplicity of the evolution in (^ is that the eigenfunctions corresponding to the states 
with lowest energy coincide with the coefficient functions $Ar entering the expansion in 
Eq. (0) m m 



$Ar(x,) = vI/W^(x,), q^ = {N + 3Y + 3/8, q^'^ = {N + 3 T 2)^ - 53/8 . (21) 

To be precise, (xi) coincides with the ground state of H^^l and ^^'^(xj) with the two 
lowest eigenstates of Ti,!^^ with opposite parity with respect to permutations of quarks, 
Xi ^ X3. The corresponding eigenvalues determine the lowest anomalous dimensions for 
each A^, mmqSM^g = En^o — ^cEn + 0(1/ Nc) and they are equal to 



E^'^ = 2^(Ar + 3) + l^-l + 27, 



E ' 



E'^ = 2^{N + 3) + - ^ + 27e • (22) 
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We will often omit the subscript "g = go" for these special states. Note useful relations: 

2E'^ = El'+ + El'- , 2$^ (x,) = <l>^'+(x,)±$^'-(x,), ($^-|$^+) = 0. (23) 

Since the wave functions of eigenstates with different energies are mutually orthogonal, 
Eq. (^) implies that the single term q = qo survives in the sum (^) in the leading large 
Nc limit. As a consequence, the moments of twist-3 parton distributions are entirely given 
by the reduced matrix elements of the corresponding local operators Oiy^g^, Eq. (|13[) P, 0]: 

■1 

(iV + 4) / 'dxx'^^^hLix.Q^) = {{O^-m)) , (24) 



N + 3 



N + 2J_, 

Combined with the scale dependence of matrix elements, Eq. (|14[), these relations are 
equivalent to the DGLAP evolution equations in (|) with NcEj^ = — dz z^^'^pj^\z). 

Treating the ~ 1/Nc contribution to the Hamiltonian in (|]) as a perturbation, we 
expand 

SN,q = N,EN,g + K^5EN,q + . . . , (25) 
^iv,g = + N:^5^>N,q + . . . , 
with the usual quantum-mechanical expressions 

5E^,, = -2||vI/(;)j|-2(xl/(;)j7^«|v^W^^ (26) 

/^(O) 1-17(1)1^(0) . ^(0) / X 

To this accuracy, moments of the twist-3 distributions are no longer renormalized mul- 
tiplicatively and have to be expanded in contributions of multiplicatively renormalizable 
operators, Eq. (|T^: 



Cf{N) j dxx^^^ fix, Q-) = {{0^,UQ'))) - ^ E 



(28) 

where Ce = 1, c,,^ = + c^^ = (iV + 3)/(A^ + 2), <l>iv = ^JJ^^^^^ is one of the functions 

$jv'~- ^N,qoiQ^) is the quark-gluon-antiquark operator with the lowest anomalous 
dimension corresponding to the wave function '^N,qo = ^7v,go ~'~^c~^'^^Af,i?o ^^'^ normalized 
at the scale Q^. Thus, the calculation of corrections to the moments of the structure 
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functions consists in two separate tasks. First, one has to calculate the 0{1/Nl) cor- 
rections to the anomalous dimension of OAr,go) or equivalently to the energies (|22D of the 



lowest eigenstates and, second, calculate the mixing coefficients 



(0) 
N,q 



\n^M^j with 



higher levels; evolution of the latter can be taken into account in the leading \axge-Nc 
approximation of Refs. [I 



np. We address both questions in what follows. 



4. The 0{1/Nl) correction 5En to the energy of the lowest eigenstates is given by 
Eq. (^) for q = qo- The calculation is most easily done by going over to the so-called 
conformal basis p, ^, [1^, |1T|. The idea is to define a basis of (orthogonal with respect 
to (p!5D) polynomials that satisfy the conformal constraints in ([T^) and, in addition, 
diagonalize the two-particle Casimir operator (H) in one of the channels, e.g. 



LuYnu^ {xi,X2,X3) = (n 5/2)(n + 3/2) Y^^^^^(xi,X2,X3) 



(29) 



and, similarly, Yj^^'^^^ 



and yj;^^'^. Here, the superscripts indicate the order in which 



the conformal spins of the three partons sum up to the total conformal spin N + 7/2. 
The three different sets of F-functions are related to each other through the Racah 6j- 
symbols of the SL{2) group. Permutation symmetry between the quarks implies that 
YNn^{xi,X2,X3) = (— 1)" y^^^^^(x3, X2, xi). Tlic explicit expressions for the F— functions 
are given in terms of the Jacobi polynomials p!0| . 
The expansion coefficients of the lowest eigenstates 



ISD in each of the three conformal 



basis are easily obtained JlT], |l^] 



N 



*5^ = E(-i)"(^+^+5) 



ra=0 
N 



N + 3 
n + 1 



± (n + 3) 



5^[(-l)^-"±l] (2n + 3)(n + 2^^^^ 



^(12)3 
^ N,n 



n 



(30) 



n=0 



9^ (31)2 

2(Ar-n) + 3Tl 



and = ($^' -|- $jy' )/2. The normalization is 

(12)3,12 _ in + l){N -n + 1) 



\\Y, 



N,n 



\\Y, 



(31)2 II 2 



{n + 2){n + 3){N + n + 5) ' 
2(n + l)(A^-n + 1)(A^ 



n + 2) 



7V,n 



(31) 



(n + 2){2n + 3){N + n + 4:){N + n + 5)' 

Notice that the two-particle kernels defined in (|^) become diagonal in the corresponding 
basis. Expanding each contribution in a suitable basis, one obtains the matrix elements, 
{^'j^^IH't'I'^n ) and (^^ IH^'^I^^ ), as finite sums over two-particle spins n. Full 
expressions are rather cumbersome and will be presented elsewhere ||12|| . 

Expanding the resulting expressions for in powers of 1/ ( + 3) we obtain 



-2(ln(iV + 3)+7, + ^-y l+O 



6E 



N 



N + 3 



3- 



ln2(iV + 3) 
(AT + 3)2 

5 



(32) 



(ln(iV + 3) + 7j - 



TT 

y 



o 



ln^(Ar + 3)' 

{N + 3y , 
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1 


2 


10 


20 


50 


100 


exact 
4'+: Eq.(P 
la,Tge-Nc 


6.1111 

6.1146^ 

6.5000 


8.1111 
8.1112 
8.7500 


9.5902 
9.5805 
10.4000 


15.3901 
15.3779 
16.8885 


18.6260 
18.6193 
20.5144 


23.2234 
23.2233 
25.6717 


26.8234 
26.8260 
29.7134 


exact 
large- Ai'c 




11.5556 

10.9640^ 

11.7500 


12.2111 
11.8972 
12.8000 


16.3820 
16.3258 
17.8116 


19.1948 
19.1742 
21.0362 


23.4791 
23.4763 
25.8981 


26.9596 
26.9612 
29.8299 


exact 
S^: Eq.® 
large- iVc 


8.5556 

7.9794^ 

8.5000 


9.7550 
9.5376 
10.2500 


10.8914 
10.7389 
11.6000 


15.8758 
15.8519 
17.3500 


18.9033 
18.8968 
20.7753 


23.3480 
23.3498 
25.7849 


26.8932 
26.8936 
29.7716 



Table 1: The lowest anomalous dimensions in the spectrum of twist-3 operators for 
different calculated by taking into account the 0{1/Nc) correction, Eqs. (|32D and 
(p5|), in comparison with the corresponding exact numerical results and leading large- A^^c 
expressions (^). 



With this correction, Eq. (|25D gives an excellent description of the lowest anomalous 
dimension in the spectrum of twist-3 operators for all integer > 0, see Table 1. Note 
that to this accuracy £j^'^ + £j^'~ = 2£^, cf. Eq. (p3|). | 
The following comments are in order. 

First, we note that 6En has the same large- A^ behavior, SE^ ~ InA^, as the leading 



large- A'c result in (|22|) . The coefficient in front of the InA^ term at large A^ is redefined, 
therefore, from 2A'c to iCp = 2(A'^ — l)/Nc, in agreement with |Q ^. It is possible to 
show |]T^ that this coefficient is exact and higher order 1/A^^ corrections do not grow 
with N, £N,go = '^Cf \nN + 0{N'^). 

Second, the constant 0{N^) term in the large A^ expansion of the anomalous dimen- 
sions (^) does not agree with 0, H, ^. The difference is due to the contribution of the 
operators Vql^ and u'^g in (H) that was overlooked in the previous works. The result for 
the 0{1/N) correction is new. 

Finally, from the expressions in ( p^ ) it is easy to read out the corresponding modifi- 



cation of the DGLAP splitting functions (g), Pf{z) = Pf\z) + P'f'\z). For example 



(33) 



This result is exact to the 0{1/N^) accuracy and neglecting all terms that vanish at 
z ^ 1. Expressions for Pe'^\z) and P'j^^{z) have similar structure. 

The opposite limit of the small— 2; behavior of the DGLAP splitting functions Pf{z) is 
of special interest. It is well known that this behavior is governed by singularities of the 



^The difference of these values with the exact result is entirely due to the truncation of the 1/(A^ + 
3)— expansion in (|3^). Since only one independent operator exists in these three cases, there is no mixing 
and the 0{1/Nc) approximation to the energies £^^q, SJ/^q and Sj/Zi is in fact exact. 
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anomalous dimensions 8n in the complex iV-plane. In particular, the leading large- iVc 
anomalous dimensions (p2D have simple poles at negative integer < —3. The pole at 
= —3 corresponds to P^^\z) in (Q). With this connection in mind, we have 

studied the analytic structure of corrections to the anomalous dimensions, 5E]^, 

continued analytically to the complex N plane. The calculation is straightforward, albeit 
tedious |T2[. We find that the analytic continuation has to be performed separately for 



even and odd A^, as familiar from studies of the evolution of twist-2 parton distribution 
beyond the leading order |]14[. Similarly to the latter case, we define[| 

1 1 + (-1)^ 1 - (-1)^ 

dxX^f{x) = / dx fevenix) -\ / dxX^ fodd{x) (34) 



1 ^ Jo ^ Jo 

and consider the evolution of even- and odd-signature component of f = e , ,g2 sepa- 
rately. When the restriction to even (odd) N is imposed, the normalized matrix elements 
($Ar|7i'^^^|$7v)/||$Arp dcvelop singularities at, generically, N = 0,-1,-2. With one ex- 
ception, all singularities to the right of = —3 are simple poles and appear due to 
vanishing of the norm of the leading large- A^c wave functions ||$Arp. They give rise to 
small-^ behavior of the correction to the DGLAP splitting functions of the form 

,(1) , , .-o 1 r . 17 - 2irM 



This asymptotics is more singular than that of Pg2\z) and is in apparent contradiction 
with the Regge theory expectations. The origin of these "spurious" singularities and their 
physical significance is not clear and deserves a special study that goes beyond the tasks 
of this letter. We would like to stress that appearance of singularities of the anomalous 
dimensions to the right from A^ = —3 can well be artifact of the 1/Nc expansion, since 
close to the singularities the correction dominates over leading order term in ( [25| ) 

and the expansion formally breaks down. 



Another delicate point is that analytic continuation in A^ has to be done in (|28| ) 



simultaneously with the analytic continuation of the sum over anomalous dimensions 



£N,qe belonging to different "trajectories" parameterized by integer i |]TT|. To the best of 
our knowledge, the problem of analytic continuation from a set of discrete points SN,q on 
a (A^,g)— plane has never been addressed in mathematical literature. 

To illustrate that the problem is nontrivial, consider the trace of the (full) Hamiltonian 
(I) over the subspace spanned by wave functions with given A^ > 



N 



TTNn = J2^N,q,. (36) 



e=o 



"This corresponds to the decomposition over partial waves with definite signature. 
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It is equal to the sum of all + 1 anomalous dimensions and is easily calculated exactly 
as a sum of diagonal matrix elements in any suitable basis. We obtain, in particular 



TtnH 



(0) 

s+ 



2(2iV + 5) ij{N + 3)+-f 
{N + 2) U(iV + 3) + 7, 
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2 N + 2 N + 3 

iV-5 + — ^ r + -ln2 

2(A^ + 2) 2 



(37) 



+ (-1 



4 



^ , N 3 
^•t + 2 



+ 



1 



+ 



2(A^ + 2) 2(A^ + 3) 



Notice now that the trace of already the leading large-A^c Hamiltonian Ti.^'^^ is singular 
at iV = —2. This implies that either one (or several) trajectories of the anomalous 
dimensions, SN,qi, becomes singular at this point, or the (analytically continued) sum 
over the trajectories in ( pSj ) diverges. In both cases, analytic continuation of the lowest 
trajectory ( ^2]) beyond = — 2 is questionable and, therefore, asymptotic expressions 
in ( p5| ) (or, equivalently, the DGLAP-evolution for z < have to be taken with 

caution. 



According to (^ 



the moments of twist-3 distributions receive contributions 
from the whole tower of operators with non-minimal anomalous dimensions at the level 
of 0{1/N^) corrections. Aside from the "energy denominators", the coefficients in front 
of these operators in (p8D are given by (normalized) matrix elements (£'|7i^^-'|£ 
the perturbation Ti,^^^ 



0) of 



II " ' 



{3t 



N,qi, II II ^ N,qe II 

where £,£' = 0,1,..., A^ numerate the anomalous dimensions from below. Our main 
observation is that the matrix elements ( PH] ) are concentrated close to the diagonal i = £', 
see Figs. 1, 2 and Table 2. In addition, the lowest anomalous dimensions E^^q^ = En 
(p^) are separated from the rest of the spectrum by a finite gap E^^g^ — E^ = 0.227 as 
A^ — s> oo p, |lT|. These two properties guarantee that the expansion in (pSD is rapidly 
converging and remains well defined in the large- A^ limit (when the number of eigenstates 
diverges). The general structure for the cases of 1-Ls+ and TIt is very similar; we present 
the results for as related to the structure function g2{x) and therefore being of more 
direct phenomenological significance. 

At large A^ the matrix elements (^IT^'^^-'IO) can be calculated semi-classically using the 
approach developed in 0]. It turns out that the matrix elements |(^|7i'^^^|0)| approach 
maximal value ~ 1/Vln A^ at £ ~ InA^, while for £^\n.N they rapidly decrease as l/£^. 
As the result, the sum in ( PBD is effectively cut off at £ ~ In A^ terms. For example, for 
A^ ~ 100 only ~ 5 first (lowest) levels give a sizeable contribution to the r.h.s. of (pSf), 
cf. Fig. 2. Analytic expressions for the matrix elements {£\H^^^Q) in the large A^ limit 



are complicated and will be given elsewhere |12 . 
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20 40 60 80 100 

Figure 1: The density plot of the mixing matrix \{e\nf+\€)\ &t N = 100; £ = corre- 
sponds to the lowest eigenstate of the large- A'^c Hamiltonian H^"!- Notice "chess- 
board" structure with alternating large (dark) and small (light) elements. 

An independent argument in favor of the dominance of the operators with lowest 
anomalous dimensions comes from the structure of matrix elements {{ON,q)) in the large- 
N limit. In an analogy with twist-2 distributions, we can write the matrix elements 
{{ON,q)) as (generalized) moments of nonperturbative (chiral-odd or chiral-even) distribu- 
tion function Df{xi,X2, x^) describing quark-gluon-antiquark correlations in the nucleon, 
schematically^ 

{{ON,q)) = j [dx]'^N,q{Xl,X2,X3) Df{xi,X2,X3), (39) 

where [dx] = dxi dx2 dx^ 5{xi + X2 + X3) and the variables Xj have the meaning of the 
quark, antiquark and gluon momentum fractions, — 1 < Xj < 1. 

At large A^, the eigenf unctions '^N,q{xi) are sharply peaked at the boundary of the 
integration region xi + X2 + = so that the matrix elements (^) probe the behavior 
of the quark-gluon distribution functions D{xi) near the kinematical boundaries. The 
precise position of the peak depends on q. It turns out that the eigenf unctions of the low- 
lying states are peaked at xi = — X3 = ±1, X2 = corresponding to configurations when 
the quark and the antiquark carry all the momentum of nucleon and the gluon is soft. 
For higher states the position of the peak is shifted gradually to 2xi = 2x3 = —X2 = ±1 
corresponding to a hard gluon and (relatively) soft quark and antiquark. If one assumes 
that the nonperturbative distributions are generated by gluon radiation already at low 

**see Appendix A in ||ri| for exact expressions. 
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Fieure 2: Exact numerical results for the coefficients v-f — '"'^im,^. „ that enter the expansion 

(p!7|) of the = 100th moment of ^'2(2;, Q^) in contributions of multiplicatively renormal- 
ized operators. Notice smallness of all coefficients with ^ > 1 and their different behavior 
for even and odd 



scales, the configurations with soft gluons should be enhanced and those with hard glu- 
ons suppressed. We conclude that the contribution of the conformal operators with large 
anomalous dimensions to the moments of the structure functions, ([2^), is suppressed 
because of smallness of both perturbative mixing coefficients and the corresponding non- 
perturbative matrix elements as describing rare parton configurations in nucleon involving 
a hard gluon. 

To summarize, we have argued that (up to 0{1/N^) corrections) the contribution 
of multiplicatively renormalizable operators with higher anomalous dimensions to the 
moments of twist-3 structure functions e(x), hi and g2{x) is small as compared with 
that of the ground state operator. With the increasing precision of the experimental 
data this contribution can be estimated in a "two-channel" approximation in which one 
supplements the dominant ground state quark-gluon component with one extra effective 
partonic component in order to account for the admixture of ~ InA^ lowest eigenstates. 
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